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We analyze the near-field interaction between the resonant sub-wavelength elements of a meta- 
material, and present a method to calculate the electric and magnetic interaction coefficients. We 
show that by adjusting the relative configuration of the neighboring split ring resonators it becomes 
possible to manipulate this near-field interaction, and thus tune the response of metamaterials. 
We use the results of this analysis to explain the experimentally observed tuning of microwave 
metamaterials. 



I. INTRODUCTION 

Metamaterials, which are typically regular arrays of 
sub- wavelength resonant particles, offer us a new degree 
of freedom in controlling the electromagnetic response of 
matter. Thus we are no longer completely constrained 
by the properties of existing materials, but can tailor 
the response in an almost arbitrary fashion, for exam- 
ple achieving very higlfD very lo\\0, and negative^ val- 
ues of refractive index, permittivity and/or permeability. 
Because of the inherently strong dispersion of resonant 
metamaterials, they must be modified in order to oper- 
ate in a different frequency band. Therefore, there is a 
significant push to have a further degree of control over 
these materials — tunability of their response. 

Fortunately, the engineered nature of metamaterials 
allows their properties to be controlled externally, either 
by dynamically modifying their structure, or by adding 
some nonlinear inclusion and controlling with external 
fields^. Examples of the latter approach include the in- 
troduction of varactor diodes^, ferroelectric^ and photo- 
conductive semi-conductors^. On the other hand, even 
without resorting to such exotic (and often lossy) con- 
stituents, there is a great deal of freedom to manipulate 
the structure itself, and this is the approach we take here. 
We consider specifically the split ring resonator (SRR) as 
one of the most important metamaterial elements, not- 
ing that whilst the fine details of near-field interaction 
are structurally specific, our approach can be applied to 
a wide variety of structures. 

An analytical model for the magnetic response of a 
sub-wavelength array of identically oriented wire loops 
loaded with a capacitance^ takes into account the mu- 
tual interaction of all the elements in the lattice, which is 
essential for deriving the effective permeability correctly. 
Although that analysis is limited to the quasi-static case 
accounting only for magnetic near-field interactions, it is 
crucial for revealing the consequence of lattice changes. 
These tend to be overlooked by otherwise rigorous ap- 
proaches which include spatial dispersion but develop 
Lorentz local field approaches, based on nearest-neighbor 
interaction^ or point-dipole approximation^. 




FIG. 1: (Color online) Several approaches to modify the lat- 
tice for metamaterial tunability: (a) a change of the lattice 
constanlp^, (b) a continuous shift of the layer and (c) a 
super-lattice of alternating shifts of layers^. 



In particular, it was pointed out in Ref. [5]that the res- 
onant frequency of the metamaterial permeability can be 
altered by varying the lattice constants without chang- 
ing the structural units. This scheme is illustrated in 
Fig.JlJa), and has been verified by experiments in Ref.fTTl 
However, a practical consequence of this change in lat- 
tice constant is that the sample size also changes corre- 
spondingly. More recently, an alternative approach was 
suggested in Ref. [T3J introducing a shift between layers 
in order to create a monoclinic lattice, with the shift in- 
creasing linearly between layers, as shown in Fig. [ijb). 
This configuration keeps the density of elements within 
the metamaterial constant, and can tune the coupling 
between neighboring particles to modify the response of 
the complete metamaterial. However, for finite size sam- 
ples, this shift inevitably results in a significant change 
in sample shape. Thus, for practical purposes, we have 
proposed a super-lattice type of geometry, whereby only 
every second layer is shifted by the same amount, as 
shown in Fig. |ljc). This tuning scheme proved to be 
robust, and allows significant manipulation of the reso- 
nant frequency with only a small change in the sample 
geometry^. Therefore, the sample geometry and its ef- 
fective properties can be engineered almost independently 
to achieve the desired manipulation of electromagnetic 
waves. 

However, as we demonstrate below, this structural tun- 
ing of metamaterials depends very strongly on the nature 
of the near-field interactions. Since metamaterial ele- 
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ments such as split ring resonators are usually not highly 
symmetric, the relative orientation of particles within the 
lattice is of key importance. This effect is not described 
by existing circuit theory models, and can give rise to 
some surprising experimental results, which we present 
here. 

In order to understand the coupling mechanisms and 
how they are affected by the lattice shift, it is useful 
to start with the simplest geometry — a pair of split 
ring resonators. Several authors have conducted numer- 
ical and experimental investigations of coupling between 
metamaterial elements for different types of ele ments and 
relative orientation between them, in microwav e 1 13 ^ 14 ! and 
optical frequency ranges, (see the overview in Ref . IT51 and 
references therein, as well as Refs. [TBHUl) A detailed 
study has previously been undertaken on tailoring the ge- 
ometric arrangement of a pair of coupled one-dimensional 
SRR arrays to engineer the dispersion curves of magneto- 
inductive wavea^SJ 

As we have shown recently-^, a rough qualitative un- 
derstanding of the structural tuning can be achieved by 
using circuit theory with purely inductive coupling be- 
tween SRRs. However, this approximation fails to pro- 
vide a quantitative understanding of the metamaterial 
tuning in question and, most importantly, does not ex- 
plain the observed strong influence of the relative SRR 
orientation. Therefore, below we develop a new model 
based on the calculation of the fundamental mode of 
a single resonator. The knowledge of the current and 
charge distributions within the mode allows calculating 
the coupling constants of a pair of split rings. These 
constants are then used to explain the experimentally 
observed tuning response of our metamaterial samples. 

In Section |ITJ we develop our approach to calculate the 
metamaterial coupling, including a discussion of the lim- 
itations of other purely analytical methods. In Section 
III[ we apply our approach to the study of interaction 
between a pair of split ring resonators which are shifted 
laterally relative to each other, and explain quantitatively 
how the shift affects the position of the fundamental res- 
onance. In Section |IV| we apply these results to a bulk 
metamaterial and identify the mechanisms at work in the 
experimental tuning of a metamaterial slab in a waveg- 
uide. Finally Section[V| concludes the paper with further 
discussions and outlook. 



II. NEAR-FIELD INTERACTION IN 
METAMATERIALS 

Considering a single SRR, it is knowifl^ that it pos- 
sesses a discrete set of eigenmodes (standing waves) with 
corresponding cigenfrequencies. In an arbitrarily excited 
SRR, the currents and charges can be represented as a su- 
perposition of the eigenmodes. The fundamental mode 
with the lowest frequency is relevant for the magnetic 
resonance in SRRs. On the frequency scale this mode is 
well isolated from the higher-order modes, and we can 
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FIG. 2: Numerically calculated (a) current and (b) charge 
distribution across an SRR at resonance 



restrict ourselves to the single-mode approximation ne- 
glecting the excitation of higher modes. 

The time dependent charge density p and current den- 
sity J in a resonant element with excited fundamental 
mode can be written in the most general form of a stand- 
ing wave: 



p(x,i) = Q(t)g(x), 



J(x,t)=/(i)j(x), 



(1) 



(2) 



where q and J describe the charge and current distribu- 
tion in space. In SRRs, the variation of the current dis- 
tribution across the width of the conductive track could 
be neglected, however for generality our approach takes 
into account the complete 3-dimensional surface-current 
distribution. 

To satisfy the conservation of charge 



V ■ J 



dp(x, t) 
dt 



we imply that 



lit) = Q(t), 



V 



JW 



-«(x). 



(3) 



(4) 



(5) 



Thus if the the current is known, it is easy to find the 
charge distribution and vice-versa. The mode profile ob- 
tained numerically for our SRR geometry is shown in 
Fig. [2j The current j is symmetric and reaches its max- 
imum at the point opposite to the gap. In accordance 
with Eq. |5]), the charge distribution q(x) is antisymmet- 
ric and goes through zero where j(x) is maximal. We see 
that q(x) reaches its maximum magnitude near the gap. 

In the single mode approximation, the dynamics of the 
SRR can be fully described by the time-dependent am- 
plitude Q(t), and we may write the SRR Lagrangian as 
a sum of terms quadratic in Q and Q: 



C = AQ 2 - BQ Z , 



(6) 
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where A and B are constants which will be discussed 
below. Accordingly, the SRR energy reads: 



E = Q 



dC 
dQ 



C = AI 2 + BQ 2 , 



(7) 



and is nicely separated into inductive (magnetic) and ca- 
pacitive (electric) parts. Clearly, for a passive SRR we 
require A > 0, B ^ 0. 

The Lagrangian equation of motion 



d_ dC 
didQ 



dC 
dQ 



(8) 



yields that the dynamics of a single SRR is described by 
the oscillator equation for the charge amplitude: 



Q(t) + ulQ{t) = 0, 



(9) 



and the fundamental mode resonance occurs at the fre- 
quency wo = \JB I A. 

Note that in contrast to the known modifications of 
the Lagrangian formalism for mctamaterials (see e.g. 
Ref. I23|) . here we do not rely on an equivalent circuit 
model. Although one might identify the parameters 2A 
and 1 /2B as effective inductance and capacitance respec- 
tively, below we evaluate them explicitly from the known 
fundamental mode shape. In fact, the strongly inhomo- 
geneous mode profile (see Fig. [2| suggests that it is un- 
likely that the correct values of the parameters would 
agree with those calculated from a circuit analysis. Ad- 
ditionally, to find the resonant frequency, we do not need 
to calculate A and B explicitly, only their ratio. 

For the purposes of our analysis, it is sufficient to con- 
sider the case of a pair of SRRs, and it will subsequently 
be shown that this explains all of the important features 
observed in our experiments with a bulk metamaterial. 
In this case, the Lagrangian can be written as a sum of 
the single SRR Lagrangians and coupling terms, which 
we also write as quadratic in currents and charges: 

C = A{Q\+Ql+2aQxQ 2 )-B{Ql+Ql+2pQiQ 2 ). (10) 

The parameters a and (3 are the dimensionless constants 
of magnetic and electric near-field interaction respec- 
tively. 

The corresponding Lagrangian equations of motion 
yield the system of ODEs for the time-dependent am- 
plitudes Qi 2 : 



Qi + uj 2 Q 1 = -aQ 2 - PuJqQ 2 , 



Qi + ^lQ 2 



-at 



/H 2 < 



(ii) 

(12) 



Solving these equations one finds that a pair of resonators 
exhibits two resonances: symmetric and antisymmetric. 
For the symmetric resonance, Qi = Q 2 , which yields the 
resonant frequency 



1+/3 



(13) 



while the antisymmetric mode with Q\ 
frequency 



WAS = W 



1-/3 



1 



-Q 2 has the 



(14) 



The described resonance splitting is well known in the 
theory of harmonic oscillators. Generally, bringing to- 
gether two oscillators of the same resonant frequency in- 
troduces coupling between them, which results in split- 
ting into two modes. Examples have been shown of SRR 
resonant frequency as a function of some coupling param- 
eter, e.g. mutual orientatiorPS or twist angle^, and typ- 
ically demonstrate a splitting or hybridization of modes. 

As we see, the direction and strength of the resonance 
shift are determined by the coupling constants a and (3. 
To evaluate them, we use the expression for the electro- 
magnetic energy following from the Lagrangian (|10h : 



E = A(lf + Ij + 2aI 1 I 2 ) + B{Q\ + Q 2 2 + 2pQ 1 Q 2 ). (15) 

The first group of terms gives the magnetic energy and 
the second group describes the electric energy. 

A possible route to calculate a and f3 is to approxi- 
mate the electric and magnetic response of each ring by 
a few terms of the multipole expansion. The problem 
with this approach is that it is based on the assump- 
tion that the observer (i.e. the second SRR) is at a large 
distance compared to the dimensions of the source. This 
requirement is strongly violated in our metamaterial sam- 
ples, where the separation between rings is actually much 
smaller than the outer ring diameter. This is essential for 
achieving strong tuning by lattice manipulation. 

Therefore we have chosen to calculate a and j3 numer- 
ically from the known charge and current distributions, 
q{x) and j(x), of the fundamental mode in a single SRR. 
Indeed, in the single mode approximation, the energy of 
a pair of SRRs reads 

E = QtW Etll + Q 2 2 W E , 22 + 2Q 1 Q 2 W E , 12 + 

IiW H ,n + I 2 W H , 22 + 2hI 2 W HA2 , (16) 



where the parameters: 



w E , mn = / d »x / d 3,'-^M x l 



(17) 



d 3 ^oj(x)-j(x') _ (lg) 

4-7TX-X' V ' 



The integrals can be easily evaluated once the charge and 
current distributions are known. The integrations over x 
and x' are over the same ring if m — n or over different 
rings otherwise. Accordingly, V\ is a volume containing 
only the first ring, and V 2 is a volume containing only the 
second. The singular terms at x = x 1 are handled using 
the analytical formulas given in Ref. [5^1 
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Comparing Eqs. ([15]) and ( 16 1 shows that the coupling 



parameters can be evaluated as 



W H ,12 a W E ,12 



W. 



if, 11 



w, 



ESI 



(19) 



For comparison purposes, when inductive coupling is 
the dominant interaction mechanism between the SRRs 
in a metamaterial, we are able to consider an array of 
split rings as an array of current loops with some mutual 
inductance between them. This approach can then be 
used to define the effective permeability of a metamate- 
rial sample^. In particular, for thin wire loops with their 
axes oriented in the same direction, the mutual induc- 
tance can be founcPH by numerical integration: 



L nn ' (r) 



4-7T 



2tt 2-rr 



d^idt/32 cos(y>i— <£>2)X 







[p 2 +z 2 +2rl{l-cos{ip 1 -ip 2 )) + 2pr (cos(p 2 -cos<pi)] 2 , 

where the distance vector r between the ring centers has 
been decomposed into a radial component p and axial 
component z, is the ring radius and <pi and if2 repre- 
sent the angle about each ring. We can use the calculated 
mutual inductance to define an equivalent magnetic in- 
teraction parameter 



«l — L nn >/L 



(20) 



(where L denotes the self-inductance of one element) 
which should approximate the interaction energy calcu- 
lated from (19). Asymptotically this interaction decays 



as 1/r, so in a large array the nearest neighbors provide 
the strongest contribution, but do not necessarily dom- 
inate over all others. Clearly, this interaction is highly 
anisotropic^, being positive for rings on the same axis, 
but negative for rings in the same plane. 

It is also possible to develop equivalent circuit models 
to calculate the electric interaction between rings. For a 
pair of coaxial rings, the total interaction can be reliably 
modeled as a circular parallel conductor transmission 
line, as for broadside-coupled split ring resonators^ or, 
alternatively, with an extended circuit model accounting 
for the distributed capacitance and inductance^. How- 
ever, once we introduce some offset between the rings, 
these approaches are not applicable and so we do not 
consider them here. 



III. TUNING INTERACTION BETWEEN A 
PAIR OF SPLIT RING RESONATORS 

Having developed an approach for calculating near- 
field interaction between a pair of rings, we now apply 
it to a canonical system which has the basic properties 
of our experimental arrangement. We consider a pair of 
SRRs, either identically (gap-to-gap) oriented or rotated 




FIG. 3: Geometry of the pair of rings (a) broadside-coupled 
and (b) gap-to-gap orientation 



a 1.5 







/ 1 






If 


/ 1 
' / 
/ 

/ ; 










*- v J 








* 7 s 

s * v 








S * x 

X S *• » 




/Ac 




^ ^ ^ % -®gt g 
a bc\ @gtg n n *a r \ 



-0.1 0.0 0.1 0.2 

Normalized Interaction energy 



FIG. 4: (Color online) Magnetic (a) and electric (/3) coupling 
parameters for broadside coupled (be) and gap-to-gap (gtg) 
orientation of a pair of shifted SRRs. The magnetic interac- 
tion calculated from the mutual inductance is given as ax. 



by 180° with respect to each other (as in the broadside- 
coupled SRR), and subject to a lateral offset 5a. The 
geometry and incident polarization are shown in Fig. [3] 
The rings have average radius tq = 2.25mm, track width 
of 0.5mm, metal thickness of 0.03mm, gap width of 1mm, 
and are separated in the transverse direction by 1.5mm. 
The resulting resonant frequency is 10.6GHz 



We plot the interaction energy calculated from Eq. ( 19 ) 
in Fig. [4j for different offsets between the rings. It 
can be seen that the electric coupling parameter /3 is 
nearly symmetric between the two configurations. This 
can be understood from Fig. [2j where we see that the 
charge distribution has a strong dipole component ori- 
ented in the x direction. For the symmetric mode of the 
broadside-coupled orientation the charges accumulated 
on the closest sides of the SRRs have opposite signs, 
the total charge distribution has the nature of a pair 
of anti-parallel dipoles, thus We, 12 < 0, and (5 is also 
negative. In contrast, in the gap-to-gap orientation, the 
closest charges are of the same sign, the total charge dis- 
tribution becomes like a pair of parallel dipoles and the 
parameter (3 > 0. 

At an offset of about one ring radius (5a « rn), the 
charge on one ring is approximately equidistant from 
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both positive and negative charges on the opposite ring, 
thus the net coupling passes through zero. At larger off- 
sets, the electric coupling changes its sign but remains 
smaller, since only the nearest halves of the SRRs are 
effectively interacting, with this interaction decaying to- 
wards zero as the offset increases. 

The magnetic interaction energy is also quite different 
for the two orientations, with the magnetic interaction 
Qi calculated by Eq. ( 20 ) lying between the gap-to-gap 
and broadside curves. For the broadside-coupled case, 
the situation is qualitatively similar to the analytical re- 
sult «i. At low offset, the magnetic field of one ring 
cuts through the other ring in the same direction to the 
surface normal, thus reinforcing the magnetic field and 
increasing the total energy. As the offset is increased, the 
situation gradually shifts to become like a pair of loops 
in the same plane, where the field from one ring cuts 
through the other in the opposite direction with respect 
to the surface normal. Hence, ctbc undergoes a change in 
sign. However, in comparison to the ring with uniform 
current, the coupling is substantially more negative. This 
is due to the current maxima being on opposite sides of 
the rings, and hence further away from each other. 

For the gap-to-gap orientation the magnetic interac- 
tion is much stronger for low 5a. This is due to the cur- 
rent maxima being located near each other which pro- 
duces a stronger contribution to the integral in Wh,12 
thus increasing a. As the rings are further separated from 
each other, the interaction energy reduces, but does not 
undergo a change of sign. We can intuitively understand 
this by neglecting the small contributions of the current 
in the region near the gaps, thus we effectively have two 
linear current elements in the same plane which always 
interact with the same sign. However, this balance is not 
universal and is determined by the specific geometry and 
parameters. To check this, we studied a geometry with a 
very small gap, so that the current distribution was much 
more homogeneous, with lower resonance frequency. In 
this simulation (not shown) the magnetic coupling did 
change sign and both values of a converged closely to 
a L . 

In order to verify that the calculated coupling correctly 
describes the frequency splitting of this system, we com- 
pare the frequency shift predicted by Eq. (13) and (14) 



with that obtained from the full numeric simulations. For 
consistency with our interaction energy approach we as- 
sume a homogeneous free-space background. We use the 
frequency domain solver of the commercial software pack- 
age CST Microwave Studio®! to model a pair of rings, in 
a unit cell with periodic boundary conditions in the di- 
rections transverse to the propagation direction. This 
periodic system enables us to define a transmission co- 
efficient, and the boundaries are 10mm from the rings. 
This value is chosen to be large enough so that there 
is no significant interaction with periodic neighbors, yet 
small enough to avoid significant scattering into higher 
order diffraction modes. Thus we can consider this the 
limiting case of a highly dilute metamaterial slab. 
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FIG. 5: (Color online) Numerical results. Transmission spec- 
trum for a pair of (a) broadside-coupl ed a nd (b) gap-to-gap 
oriented rings. Solid line: ujs from Eq. (131, broken line: lj as 
from Eq. (Til) 



The transmission spectrum as a function of offset is 
plotted in Fig. [5] We see that all the important fea- 
tures of the mode splitting are represented correctly by 
our single-mode theory of coupled SRRs. Expansion of 
Eqs. (13) and (14) for small coupling predicts frequency 
splitting of Aw = ±(/3 — a)/2, hence the curves are ap- 
proximately symmetric about ujq. 

The strong splitting observed for the broadside- 
coupled orientation is due to the opposite signs of the 
electric and magnetic coupling. At Sa/r w 1.1 both a^ c 
and ftbc change signs, hence the crossing of the symmetric 
and anti-symmetric modes is observed. In contrast, we 
see that for the gap-to-gap orientation, for small offset 
Oig t g and f3 g tg are of the same sign, and thus they have 
an opposing effect, resulting in small frequency splitting. 
Since j3 decreases much faster and changes sign, the result 
in maximum frequency splitting for 5a /r$ between 1.1 
and 1.5. We note that in Ref. Q2]the resonant frequency 
based onax was compared with experimental results for 
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the gap-to-gap orientation, and strong disagreement was 
found. 

It can be seen that the calculated transmission through 
the cell exhibits different depths of the resonance for 
the symmetric and antisymmetric modes. This is due 
to the different efficiency of coupling between the modes 
and the incident plane wave. For instance, it is impossi- 
ble to excite the anti-symmetric mode with a normally- 
incident plane wave for 5a = 0, since both rings are ex- 
cited in-phase. As the offset is increased, some retar- 
dation between the rings occurs, and excitation of the 
anti-symmetric mode is allowed. 

There are several reasons for the small quantitative 
disagreement between the exact calculations and those 
based on the calculated interaction energy. Firstly, the 
minimum of transmission occurs at a frequency slightly 
different from the resonant frequency, due to coupling ef- 
fects (impedance matching) between the incident wave 
and the ring. Secondly, there may be some small con- 
tribution of higher SRR eigenmodes due to perturbation 
of the charge and current distributions. Thirdly, there 
may still be some small influence of the periodic bound- 
aries. Finally, our developed relations neglect retarda- 
tion, which is strictly valid in the sub- wavelength limit, 
whereas the the outer radius of the rings is 0.18A at u>o- 
Retardation has previously been shown to modify the 
interaction between SRRs through its influence on the 
dispersion of magneto-inductive waves in arrays^H^U. 

We emphasize that our approach developed in Section 
[IT| is advantageous over the direct numerical calculation. 
Firstly, once the mode profile is known, calculation of the 
frequencies in Fig. [5] takes approximately 30 seconds on 
a single CPU, whereas the direct calculation of the full 
spectrum takes several hours on a multi-core machine. 
Secondly, we are clearly able to demonstrate the nature 
of the coupling, which yields insight into the tuning be- 
havior. 



IV. TUNING INTERACTION IN A BULK 
METAMATERIAL 

We now apply our approach to explain experimental 
results for tuning the response of a slab of metamaterial. 
The metamaterial is fabricated using photolithography 
to etch copper tracks onto FR4 printed circuit board, 
using the same geometry as in our numerical simulation 
of a pair of rings. The fabricated sample has 30 lay- 
ers, each with 5 rings in the propagation direction and 
is only one ring in height (i.e. a 5 x 30 x I array). The 
longitudinal period of the sample is 7mm, the transverse 
period is dictated by the sample thickness (1.5mm) as 
there is no spacing between boards. As with the pair 
of rings, we have assembled slabs with two relative ori- 
entations of the split rings in adjacent planes — gap- 
to-gap and broadside-coupled. The sample is placed in 
the center of a WR229 rectangular metallic waveguide, 
with dimensions 58. 17x 29.08mm, excited at its dominant 



(a) Broadside-coupled orientation 




(b) Gap-to-gap orientation 




f (GHz) 



FIG. 6: (Color online) Experimental transmission while tun- 
ing 8a for split ring resonator slab in waveguide, for (a) 
broadside-coupled and (b) gap-to-gap orientation of adjacent 
layers. 



TE 10 mode. We removed the influence of the coaxial 
adapters and feeding waveguide sections by performing 
a TRL (through-reflect-line) calibration^. In Fig. [fj] we 
show the experimentally obtained transmission through 
each slab of closely coupled SRRs, with the correspond- 
ing reflection shown in Fig. [7j 

The shift of resonant frequency shows good qualita- 
tive agreement with the results for a pair of rings pre- 
sented in Section |III[ with very similar changes in the 
spectrum observed (noting that Sa/a = 0.5 corresponds 
to Sa/r = 1.56). However, for the gap-to-gap orien- 
tation, numerical simulation of a system of two boards 
with 5 rings each, and periodic boundaries in both trans- 
verse directions (not shown), is qualitatively similar to 
the experimentally observed results but quantitatively 
highly inaccurate. The reason turns out to be the loss of 
symmetry when the system is placed inside the waveg- 
uide, because the upper and lower waveguide walls do not 
correspond to periodic boundaries, but instead represent 
planes of mirror reflection. Therefore, this system must 
be described as having a super-lattice arrangement in the 
vertical as well as horizontal planes, with each super-cell 
consisting of four SRRs. This cell has alternating orien- 
tation of the SRRs in the vertical direction corresponding 
to the planes of mirror symmetry, as shown in Fig. [8](b) . 
Once this unit cell is taken into account, numerical sim- 
ulations are in a good agreement with the experiment 
(Fig- [lb)). 

Naturally, numerical simulations for the broadside- 
coupled orientation also agree well with the experiment 



7 



a) Broadside-coupled orientation 




(a) Broadside coupled orientation 



3.5 4.0 4.5 

(b) Gap-to-gap orientation 




FIG. 7: (Color online) Experimental reflection while tuning 
5a for split ring resonator slab in waveguide, for (a) broadside- 
coupled and (b) gap-to-gap orientation of adjacent layers. 




FIG. 8: Schematic of the effective super-lattice geometry cor- 
responding to the waveguide measurement for (a) broadside- 
coupled and (b) gap-to-gap orientation. Dashed lines show 
planes of reflection symmetry, and the shaded region shows 
the super-cell with 4 SRRs. 



(Fig.^a)). In this case a very similar result is provided 
with simple periodic boundary conditions (not shown). 
For this orientation the super-lattice effectively formed 
by the waveguide shown in Fig.[8]ja) does not have an es- 
sentially different symmetry to the original super-lattice. 

We can conclude that in both cases the dominant mode 
of the slab corresponds to the dominant symmetric mode 
of a pair of rings, with a similar pattern of resonant fre- 
quency vs. offset occurring. The weaker coupling to the 
modes for the gap-to-gap orientation is due to the shape 
of the symmetric mode. Its magnetic field has a large 
component parallel to that of the waveguide mode, how- 
ever its electric field is primarily longitudinal, in contrast 
to the transverse electric field of the waveguide mode. 

In Fig. |9jb), we see two higher order modes, which 
most likely correspond to the higher order resonances ob- 




FIG. 9: (Color online) Numerically obtained transmis- 
sion spectrum of metamaterial in waveguide, (a) broadside- 
coupled and (b) gap-to-gap orientation. The white line indi- 
cates experimentally obtained resonant frequencies. 



served in the experimental results in Fig.[6jjb). From the 
numerical simulations we observe that the current distri- 
butions of these modes are symmetric, thus they corre- 
spond to higher-order modes of the metamaterial slab, 
and not to the anti-symmetric mode of a pair of rings. 
In contrast, in Fig. [9^a) there is a weakly coupled anti- 
symmetric mode, which we verified by inspection of the 
currents. This mode also qualitatively agrees with the 
corresponding mode of the pair of rings, with somewhat 
weaker coupling due to the increased mismatch to the 
incident waveguide mode. This mode may correspond 
to some of the smaller features observable in Fig. |6ja), 
however due to the size of these features this cannot be 
reliably determined. 

We do not consider offsets greater than 0.5a, since in an 
infinite lattice only shifts between and 0.5 are unique, 
while in a finite structure, larger shifts result in very ir- 
regular boundaries. Note that in the simulations we have 
neglected the effect of the mode profile of the rectangular 
waveguide, which would correspond to an effective vari- 



ation of the angle of incidence of the plane wave as a 
function of frequency, which can result in a different re- 
sponse due to the anisotropy and non-negligible spatial 
dispersion of the medium^. 

Clearly the coupling in the complete lattice is much 
more complicated than in the simple two-ring system, as 
the interactions between a large number of rings must 
be taken into account. In principle it is possible to ex- 



tend the analysis of Section III to an arbitrary number 



of rings. However the qualitative agreement between the 
experimental results and the modelled pair of rings sug- 
gests that the phenomenology developed for the two rings 
is generally applicable and leads to correct predictions. 

Clearly, the coupling in the metamaterial lattice is 
much more complicated than in the simple two-ring sys- 
tem, as the interactions between a large number of rings 
must be taken into account. The qualitative agreement 
between the experimental results and the modelled pair 
of rings suggests that the phenomenology developed in 



Section III for the two rings can be extended to an arbi- 
trary number of rings. 

Although an accurate generalization of our modelling 
approach to a bulk system lies beyond the scope of this 
paper, it is clear that the resulting homogenized effec- 
tive metamaterial parameters will exhibit similar tuning 
pattern due to the resonance shift. Note that the in- 
troduction of the effective parameters is justified when 
the ratio of the unit cell size to the incident wavelength 
is small. Therefore, when considering modifications of 
the lattice which create a super-lattice structure, the 
size of the super-cell should be smaller than the wave- 
length. There are homogenisation approaches in the lit- 
erature (e.g. Ref. 135(1 which include unknown parame- 
ters for interaction between resonant elements and our 
semi-analytical approach would make an ideal tool for 
evaluating these constants. 



V. CONCLUSION 

We have analyzed the near-field coupling within meta- 
materials, considering both the relative orientation and 



the offset between the centers of two neighboring res- 
onators. Using a pair of split ring resonators as a simple 
model, we have shown the coupling mechanisms at work 
in our recently proposed tuning scheme, based on the di- 
rect calculation of the interaction energy. We have con- 
firmed that these mechanisms can predict qualitatively 
the performance of a realistic metamaterial structure. 
This paves a road towards a reliable design and develop- 
ment of tunable metamaterials for various applications. 

We note that the specific geometry of the split rings 
can have a very significant influence on the qualitative na- 
ture of the coupling, including cases which run counter to 
our intuitive understanding of current loops interacting 
magnetically. 

Finally, we point out that the approach developed here 
for modeling near-field effects is particularly promising 
for metamaterials scaled down to operate at optical fre- 
quencies. In the visible, the paradigm of ideally con- 
ducting metal fails and the area of applicability of circuit 
models is rather limited. In contrast, the consideration in 
terms of excitation and interaction of plasmonic standing 
waves will provide a clear physical picture. 
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